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O , Abstract 

lO ' The Centaurus cluster (z = 0.0104) was observed with the X-ray Imaging Spectrometer (XIS) onboard 

0^ ' the Suzaku X-ray satelhte in three pointings, one centered on the cluster core and the other two offset by 

, ±8' in declination. To search for possible bulk motions of the intracluster medium, the central energy of 

Hc-likc Fc-K line (at a rest-frame energy of 6.7 kcV) was examined to look for a positional dependence. 
Over spatial scales of 50 kpc to 140 kpc around the cluster core, the central line energy was found to be 
• constant within the calibration error of 15 eV. The 90% upper limit on the line-of-sight velocity difference 

is |Ai;| < 1400 kms ^, giving a tighter constraint than previous measurements. The significant velocity 
Q ■ gradients inferred from a previous Chandra study were not detected between two pairs of rectangular 

|H I regions near the cluster core. These results suggest that the bulk velocity does not largely exceed the 

C/3 . thermal velocity of the gas in the central region of the Centaurus cluster. The mean redshift of the 

' intracluster medium was determined to be 0.0097, in agreement with the optical redshift of the cluster 

^ , within the calibration uncertainty. Implications of the present results for the estimation of the cluster mass 

are briefly discussed. 
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1. Introduction its thermal velocity, the associated non-thermal pressure 

would endanger the assumption of hydrostatic ICM equi- 

Clusters of galaxies are the most massive gravitation- librium, which is usually assumed in deriving the total 

ally bound structures in the Universe, and are thought to gravitating mass in a cluster. For example, the factor 

grow into larger systems through complex interactions be- 2-3 discrepancy between the X-ray and lensing mass es- 

tween smaller systems. Signatures of such merging events timations in some objects (e.g., Ota et al. 2004; Hattori 

may manifest themselves in non-Gaussian velocity dis- et al. 1999) could be due to this effect. Therefore, the 

tributions of member galaxies, temperature and density presence/absence of ICM bulk motion has a great impact 

inhomogeneities in the intracluster medium (ICM), and on the cluster studies, including their mass estimates and 

bulk motions of the ICM. Indeed, numerical simulations cosmological applications. 

predict the existence of bulk ICM flows with a substan- We are thus encouraged to observationally constrain 

tial fraction of the virial velocity (^ 1000 kms^^ for rich the ICM motion. Such an attempt may be carried 

clusters), lasting several Gyr after each merger event (e.g., out by examining X-ray emission lines from heavy ele- 

Roettiger et al. 1996; Norman & Bryan 1999). Therefore, ments for Doppler shifts and broadenings (e.g., Dupke & 

measurements of the ICM velocity structure provide very Bregman 2001), or by utilizing kinetic Sunyaev-Zcl'dovich 

important information with which to understand the clus- effects in the radio wavebands (e.g., Sunyaev et al. 2003). 

ter dynamics. For the former method. X-ray spectroscopy of the emis- 

If the ICM has a significant bulk velocity compared to sion lines from heavy ions, particularly the 6.7 keV K-lines 
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from Helium- like iron atoms, provides the most sensitive 
way of studying such ICM motions. By observing X-ray 
spectra in the 0.5-8.5 keV range from the core regions 
of the Centaurus cluster (Abell 3526) with the Chandra 
ACIS, Dupke & Bregman (2006) recently reported on a 
large velocity gradient of 2400 ± 1000 kms~^ over a spa- 
tial scale of 100 kpc, reconfirming their previous ASCA 
measurements (Dupke & Bregman 2001). However, their 
results might still be subject to residual uncertainties aris- 
ing from temporal and intrachip gain variability, requiring 
an independent cross confirmation. Since a bulk velocity 
of 1000 kms~^ along the line of sight translates to a shift 
in the He-like Fe-K line energy of Ai? = 22 eV, we need not 
only sufficient photon statistics and good energy resolu- 
tion, but also precise instrumental gain calibrations with 
an accuracy better than ^ 20 eV (or ^ 0.3%). 

The Centaurus cluster (Abell 3526), at the redshift of 
Zci = 0.0104, is a nearby X-ray bright cluster of galaxies, 
hosting the cD galaxy NGC 4696. It has been extensively 
studied at various wavelengths. Lucey et al. (1986) found 
a bimodal velocity distribution of galaxies, corresponding 
to the main cluster (Cen 30) centered on NGC 4696, and a 
subcluster (Cen 45) centered on NGC 4709, with a veloc- 
ity difference of ~ 1500 kms~^. They hence argued that 
Cen 30 is currently accreting Cen 45. The X-ray mea- 
sured ICM temperature is ~ 3.8 keV in outer regions, and 
decreases toward the center to ^ 1.7 kcV (e.g., Ikebe et 
al. 1999; Makishima et al. 2001; Takahashi et al. 2004). 
Although the central part of this cluster is relatively re- 
laxed (e.g., Allen & Fabian 1994), at larger radii the tem- 
perature maps derived with ASCA reveal a hot spot co- 
incident with Cen 45 (Churazov et al. 1999; Furusho et 
al. 2001), supporting the idea that there is a past or on- 
going merger episode. The central region is very rich in 
metals, with a large abundance gradient (Fukazawa et 
al. 1994; Ikebe et al. 1999; Graham et al. 2006; Sanders & 
Fabian 2006). Emitting very intense Fe-K emission lines, 
the central region of the Centaurus cluster provides an 
excellent opportunity for a study of ICM bulk-motion. 

In this paper, we present an analysis of the bulk veloc- 
ity of the gas over the central 240 kpc of the Centaurus 
cluster, conducted with the X-ray Imaging Spectrometer 
(XIS; Koyama et al. 2006) on board the Suzaku satellite 
(Mitsuda et al. 2006) . The XIS has an excellent sensitivity 
at the iron K-line energies, and a very low background. In 
addition, the current accuracy of the XIS in-orbit calibra- 
tion is already high enough to enable us to make meaning- 
ful searches for possible ICM bulk motions. The Suzaku 
X-ray telescope (XRT; Scrlemitsos et al. 2006) has a spa- 
tial resolution of ~ 2' (Half Power Diameter; HPD), allow- 
ing us to study spatial velocity variations on reasonable 
scales (~a few arcminutes). We therefore test the pre- 
viously claimed large velocity gradient near the cluster 
center. 

Throughout this paper we adopt Qm — 0.3, JIa = 0.7, 
and i/o — 70 kms~^Mpc~^; 1' corresponds to 13.5 kpc at 
the cluster redshift of = 0.0104. The quoted statistical 
errors refer to 68% confidence ranges, unless otherwise 
stated. 




Fig. 1. Suzaku XISO smoothed image of the Centaurus clus- 
ter in the 0.4-10 hcV band. The three pointing images are 
superimposed, but are not corrected for the overlapping ex- 
posure. Two corners of the CCD chip illuminated by the on- 
board calibration sources are excluded. The X-ray emission is 
peaked at NGC 4696. The subcluster centered on NGC 4709 
is about 15' away from the center and outside the XIS fields 
of view. 

2. Observation 

Three pointed observations of the Centaurus cluster 
were conducted between 2005 December 27-29, in the 
Suzaku Phase-I period: one centered on the cD galaxy, 
NGC 4696, and the other two (Offset 1 and Offset2) offset 
by —8' and +8' in declination, respectively. The details 
of the observations is given in table 1 . The two onboard 
instruments, the XIS and the Hard X-ray Detector (HXD; 
Takahashi et al. 2006), were operated in their normal 
modes during the observations. The XIS consists of four 
X-ray sensitive CCD cameras: three front-illuminated 
(FI) CCDs and one back-illuminated (BI) CCD. In the 
present paper, we use data from the three Fl-chip cam- 
eras (XISO, 2, and 3), since they have higher sensitivities 
as well as lower backgrounds at the iron-K line energies, 
compared to the Bl-chip camera (XISl). As shown in fig- 
ure 1, the extended ICM emission was clearly detected, the 
center of which {a,S) = (192°. 2058, -41°. 30424) in J2000 
coordinates, is consistent within ^ 25" of the position of 
NGC 4696. 

The data reduction was performed using HEASOFT 
version 6.0.6. The XIS event lists created in the revO.7 
pipeline processing were filtered with the following crite- 
ria: the geomagnetic cut-off rigidity > 6 GV, the eleva- 
tion angle from the earth limb > 10°, and the satellite 
outside the South Atlantic Anomaly. The net exposures 
are 29.4, 33.4, and 31.6 ks for the Center, Offsetl, and 
Offset2 pointings, respectively. We selected XIS events 
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Table 1. Details of the Suzaku observations of the Centaurus cluster 



Target name 


Scquenee No. 


Date 


Coordinates* 


Roll angle 


Exposure^ 








a (deg) 5 (dog) 


(deg) 


(b) 


CENTAURUS_CLUSTER 


800014010 


2005 Dec 27 


192.2054 -41.3111 


-253.76 


29397 


CENCL.Offsetl 


800015010 


2005 Dec 28 


192.2054 -41.444 


-253.72 


33431 


CENCL_0frsct2 


800016010 


2005 Dec 29 


192.2054 -41.178 


-253.48 


31647 



Pointing coordinates in J2000. 
^ Net exposure time of the XIS sensors after data filtering. 



with grades 0, 2, 3, 4, and 6. The XIS response matrix 
files, ae_xi [0,2,3] _20060213.rmf, and the on-axis aux- 
iliary files, ae_xi [0 , 2 ,3] _xisnorm6_20060415 . arf were 
used in the spectral analysis. 

3. Data analysis and Results 

Our primary goal is to constrain the line-of-sight bulk 
velocity of the ICM using the iron-K line at 6.7 keV. For 
this purpose, the accuracy of photon energy measurements 
is crucial. We briefly review in §3.1 the current XIS cali- 
bration uncertainty, and describe in §3.2-§3.6 our analy- 
ses of the velocity structure of the ICM on small and large 
spatial scales. 

3.1. Accuracy of the energy scale 

In order to examine the accuracy of the XIS energy 
scale, we checked the following three key points: the line 
centroid energies of calibration sources illuminating some 
corners of the CCD chips, the positional gain variations 
due to Charge Transfer Inefficiency, and the agreement of 
line centroid energies among the three pointings onto the 
same sky regions of the Centaurus cluster. 

The fiducial absolute energy scale of the XIS is pro- 
vided by Mn-Ka lines from the built-in calibration sources 
(^^Fe), which illuminates two out of four corners of each 
XIS chip. For the three pointings, we fitted these lines 
with a Gaussian, and confirmed that the obtained line 
centroids of XISO, XIS2, and XIS3 all agree within 0.02- 
0.1% of the expected value of 5.8951 keV. 

Although the calibration isotope fixes the fiducial en- 
ergy scale, the gain could vary from place to place on the 
same CCD chip (intrachip gain variation), due to Charge 
Transfer Inefficiency (CTI). The CTI increases with time 
due to irradiation by charged particles, enhancing the in- 
trachip gain variations particularly along the direction of 
charge transfer (so called ACT-Y axis). The CTI char- 
acteristics of the XIS were calibrated in orbit using two 
line-rich extended sources, the Cygnus Loop and Sgr C, 
observed on 2005 November 23-24 and 2006 February 20- 
23, respectively. As a result, the gains of the four XIS 
sensors have been equalized over the CCD chips to an 
accuracy of ±0.2% (Koyama et al. 2006). These results 
have been used to produce the revO.7 data. As the three 
observations of the Centaurus cluster were carried out on 
successive days between those of the Cygnus Loop and Sgr 
C, our analysis should be little affected by any temporal 
gain variation. 



We further examined the revO.7 data of two other ex- 
tended sources, Abell 1060 and Abell 426, observed on 
2005 November 22 and 2006 February 1, respectively. The 
ACT-Y dependence of the iron line energy was then con- 
firmed to be very small, typically within ±0.2%. From 
these calibrations, we consider that the intrachip gain vari- 
ations of the revO.7 data from the four XIS sensors are 
at most ±0.2% (±15 eV) in the ACT-Y direction. The 
gain variation along the ACT-X axis is also within ±0.2%, 
as confirmed with Sgr C, AbeU 1060, and Abell 426. 
Quantitative confidence levels of the systematic error were 
estimated based on a detailed analysis of Abell 426, which 
has the most intense iron-K lines among the celestial cali- 
brators utilized here. Since the four XIS sensors are placed 
on the focal plane in different orientations with respect to 
the satellite coordinate (Koyama et al. 2006), the posi- 
tional gain variation is effectively studied by comparing 
the line energies of the same sky regions (namely different 
detector regions) on the four sensors. We thus divided the 
XIS field of view of the Abell 426 observation into 8x8 
cells of size 2'.1 x 2'.1, and estimated the la systematic 
error on the iron line energy by calculating the chi-square 

value, defined as x' = ScenSLo(-^« " {E)f I {^1 + <^lys) ■ 
Here Ei and are the iron-K line centroid energy and 
the ItT statistical error obtained with XIS-i in a certain 
cell, respectively, {E) = Sf^QE'i/4 is the mean iron-K en- 
ergy, and Sceii sums up all cells. If a systematic error of 
dsys = 10 eV (i.e., ±0.13%) is tentatively employed, we 
obtain /v = 1.1 for 144 degrees of freedom. This value 
of (Tsys is hence considered appropriate as a measure of the 
residual gain non-uniformity. 

Utilizing the overlaps between the three pointings of 
the Centaurus cluster (figure 1), we examined whether 
the central line energy of the same sky region turns out 
to be the same when measured in the different pointings 
(and hence by different regions on the CCD chip). In de- 
tail, we derived the iron- line energy of the central pointing 
data of the three FI chips from a 8'. 5 x 15' region (roughly 
coincident with southern half of the XIS field of view in 
the central pointing), and compared it with that of the 
same sky region observed in Offsetl. The results agree 
with each other within 0.1%. Similarly, the Fe-K line en- 
ergy averaged over a northern 8'. 5 x 15' region showed a 
good agreement (0.1%) between the central pointing and 
Offset2. If we divide a central 8'. 5 x 8'. 5 square region 
centered on the cD galaxy into 4x4 cells of size 2'.1 x 2'.1 
each, the energy scales are again in good agreement be- 
tween the central pointing and the two offsets; ^ 0.2% in 
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14 cells, and ^ 0.3% in 2 cells (with la statistical errors 
being ~ 20 eV). 

Based on these examinations, we decided to take 
±0.13% (±0.2%) as the systematic error on the XIS en- 
ergy scale at the 68% (90%) confidence level. The 68% 
systematic error corresponds to ±10 eV in terms of the 
iron-K line energy, or ±0.0015 in redshift, or ±470 kms^^ 
in line-of-sight velocity. 

3.2. Analysis of small-scale velocity structure 

Now that we have confirmed the XIS energy-scale is 
accurate, we proceed to the study of the velocity structure 
of the ICM in the Ccntaurus cluster on a scale of a few 
arcminutes. For this purpose, we divided the 18' x 18' 
square XIS field of view into 8x8 cells as shown in figure 2. 
Thus, each cell covers a sky area of 2'.1 x 2'.1 (=28 kpc 
X 28 kpc). Since the energy scale is consistent among the 
three XIS sensors and among the three pointings (§3.1), 
we CO- added the data from XISO, XIS2, and XIS3, and 
merged the three pointings where they overlap. Excluding 
the calibration source regions, 52 cells in total were used in 
the spectral fitting. The composite spectra were rebinned 
so that each spectral bin had more than 30 counts. 

Figure 3 shows an example of the XIS spectrum ob- 
tained in this way, displayed over the iron energy band. 
Thus, wc observe a very intense He-like Fc-Ka line, a 
weaker H-likc Fe-Ka line, and a blend between He-like 
Ni-Ka and He-hke Fe-K/3 lines. In the 5-10 keV band, 
the source count rate measured in each cell ranges from 
1.6 X 10"'^ to 2.5 X 10~^ c s^^ per XIS sensor, with a mean 
over the 52 cells being 6.9 x lO^'^ c s~^. If the energy band 
is limited to 6.5-6.9 keV, in which the He-like Fc Ka line 
emission is prominent, these count rates decrease by a 
factor of 3-5. 

Background spectra were derived by integrating events 
over the same regions on the detector planes as the 
on-source data integration, using a blank-sky observa- 
tion towards the North Ecliptic Pole conducted on 2005 
September 2-4 for 95 ks. At the cluster center, the back- 
ground surface brightness in the 5-10 keV range is only a 
few percent of that of the ICM emission, while the back- 
ground contribution increases to ^ 35 percent near the 
edge of the XIS field of view. The XIS background spec- 
trum has several instrumental fluorescent lines in the en- 
ergy band of interest, with the most prominent one being 
Ni-Ka line at 7.48 keV. However, as shown in figure 3, 
the background spectrum can be regarded as featureless 
in the energy range where the He-like Fe-Ka line from 
the Centaurus cluster is present. We subtracted the back- 
grounds from the on-source spectra, and confirmed that 
this process does not change the Fe-K line energy. 

3.3. Modeling the iron-K lines 

To accurately determine the iron-line energies, we fit- 
ted the background-subtracted XIS spectra with a simple 
model over the 5-10 keV energy range. The model we 
chose consists of a continuum represented by the APEC 
thermal emission model with metal abundances reset to 
zero (Z = 0), and three Gaussians for the major line com- 
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Fig. 2. An image showing the spectral accumulation regions. 
The central 8x8 cells, each with a size of 2'.1 X 2'.1, are 
designated as #0 ~ # 63, and were used for the study of 
small-scale velocity structure. The two larger rectangular re- 
gions (the dashed boxes of size 7' X 18') were used for the 
larger-scale analysis. The four regions used for comparison 
with the Chandra results (Dupke & Bregman 2006) are also 
shown as PEAK-, PEAK+ (the solid boxes), R-2, and R-l-3 
(the dashed boxes). 

ponents: the He-like Fc-Ka line at 6.70 keV, the H-like 
Fe-Ka line at 6.97 keV, and a 7.83 keV line representing 
the blend of He-hke Ni-Ka hnc (7.80 kcV) and He-hkc Fe- 
K/3 line (7.90 keV). The quoted line energies refer to their 
rest-frame values. 

The He-like Fe-Ka complex, consisting of 26 lines, ex- 
hibits an asymmetric shape (e.g., Dubau & Volonte 1980). 
Although the line profile becomes essentially Gaussian 
when convolved with the XIS energy resolution [~ 160 cV 
(FWHM) at the time of the observations of the Centaurus 
cluster; Koyama et al. 2006], the Gaussian centroid en- 
ergy, denoted E'o, can no longer be identified with that 
of the strongest component, namely the resonance line 
at 6.6986 keV in the rest frame. To determine the ex- 
pected value of Eq when the ICM has no bulk motion 
and the cluster is at zero redshift. wc produced a simu- 
lated spectrum employing the APEC model and XIS en- 
ergy response. Based on ASCA measurements (Furusho 
et al. 2001; Ikebe et al. 1999), the simulation assumed an 
ICM temperature of kT ~ 3 keV, and a metallicity oi Z — 
IZq. By fitting the simulated spectrum with an APEC 
plus three Gaussian model as described above, we ob- 
tained Eq = 6.6771 ±0.0001 keV. We thus assume the rest- 
frame He-like Fe-Ka centroid energy to be Eq = 6.677 keV. 

Strictly speaking, the value of Eq is expected to depend 
upon the ICM temperature; as the temperature decreases, 
contributions from lower-ionization ions increase, making 



No. ] 

Eq shift towards lower energies. For example, the shift 
is by about 1% when kT = 1 keV. The very center of the 
Centaurus cluster is known to host a cool component with 
a temperature of kT 1.7 keV (Ikebe et al. 1999; Sanders 
& Fabian 2002; Takahashi et al. 2004). However, the 
effect of the cool plasma on the line-centroid measure- 
ment is estimated to be negligible, because the cool gas 
is localized within 1' from the cD galaxy with a rela- 
tively minor emission integral (Ikebe et al. 1999; Sanders 
& Fabian 2002; Takahashi et al. 2004), and the Fe-K 
line emissivity of the 1.7 keV gas is by more than a fac- 
tor of 2 (7) lower than that of > 2 keV (> 3 keV) gas. 
Towards higher temperatures, Eq can be regarded as con- 
stant within 0.1% as long as 2 keV <kT<8 keV. 

The intrinsic line width of the He-like Fe-Ka com- 
plex emitted from the 3 keV gas is derived to be ctq = 
27.8 ± 0.2 eV. This is smaller than the instrumental en- 
ergy resolution at 7 keV, ^ 70 eV in Gaussian sigma, 
and hence a minor contributor to the observed total line 
width. The possibility of the line broadening in excess of 
the intrinsic width is discussed in §4.3. 

3.4- Result of spectral fitting 

Using XSPEC version 11, we fitted the 5-10 keV spectra 
from individual cells with the APEC plus three Gaussian 
line model described in §3.3. The three Gaussian ener- 
gies were allowed to vary, but their relative energies were 
constrained to obey their rest-frame ratios. Since the XIS 
energy resolution 6E depends on energy E approximately 
as 6E (X E'^-^ in the energy band of interest, and varies 
gradually with time (Koyama et al. 2006), the Gaussian 
widths were included as free parameters, with their rela- 
tive values constrained to scale with the square root of the 
line energies. The normalizations of the APEC and three 
Gaussians were also left free. Although the analysis uti- 
lizes the fitting method, consistent results are derived 
if unbinned spectra are fitted by the maximum likelihood 
method. 

An example of the spectral fit is shown in figure 3. We 
obtained an acceptable fit at the 90% confidence level from 
49 out of the 52 cells (except # 2, 22, and 33), and at the 
95% level from all the cells. The reduced values, xt^ 
range from 0.68 to 1.53 in the central 6x6 cells, where 
the number of degrees of freedom ly ranges from 15 to 
132. In the outer 19 cells, xt is between 0.30 and and 
1.85, and ly between 3 and 27. The value of varies as 
the number of spectral bins changes, due to the spectral 
binning criterion. 

The line centroids have been constrained in most cases, 
with the la statistical errors of aat =3 — 4 eV for the inner 
2x2 cells, ^ 4 — 8 eV for the surrounding 12 cells, and ~ 
10 eV for outer 17 cells surrounding the central 4x4 cells. 
In most of these 33 cells, the fitted Gaussian widths are 
consistent, within the 90% statistical errors, with those 
foimd with the Mn-Ka calibration line (crinst ^ 40 eV) . On 
the other hand, those cells lying near the edge of the field 
of view have too poor photon statistics to constrain the 
line width. In analyzing the outermost 19 cells, we thus 
fixed the line width at 40 eV and the ICM temperature 
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Fig. 3. Background-subtracted XIS0-I-XIS2+XIS3 spectrum 
of cell #44 in the 5-10 keV band (the crosses), fitted with 
the APEC model with Z = for the continuum (the dashed 
line) and three gaussian lines; the dash-dot line for the He- like 
Ka, the dotted line for H-like Ka, and the dash-dot-dot line 
for a blend of He-like Ni Ka and He-like Fe K/3. The best-fit 
total model is shown with the solid line. The background 
spectrum taken from the North Ecliptic Pole observation is 
also shown (the crosses with triangle). The bottom panel 
shows the residuals of the fit in terms of the ratio of the 
residual to the uncertainty of each measurement, which is 
also the contribution to when squared. 

at 3.0 keV. 

The fitted line centroids, E'obs, are then converted 
to redshifts through z ~ {Eq — £'obs)/-E'obs using Eq = 
6.677 keV. The obtained results are plotted in figure 4, 
where the implied ICM radial velocity, Vr = c{z — Zc\) , is 
also indicated on the right-side ordinate of the panels. As 
seen from this plot, the measured redshifts are consistent, 
in almost all the 52 cells, with the optical value within 
the current 90% calibration errors (±0.002 in terms of 
redshift). 

The weighted mean (and the standard deviation) of 
the line centroid measured in the 52 cells is (-Eobs) = 
6.6133 keV (±0.0162 keV). The reduced chi-square around 
this mean is xl = J^'^^HEohs - (-Eobs))/o'fit)^ = 1-1 with 
u = 51. Thus the measured i?obs values are consistent 
with being constant at the 90% confidence level. The 
mean redshift (and the standard deviation) of the ICM is 
determined to be (zobs) = 0.0097 (±0.0024), which agrees 
with the optical redshift of the cluster Zc\ — 0.0104 within 
the calibrational uncertainty. For reference, the optical 
redshift of the cD galaxy is z^d = 0.009867 (the NED 
database). The difference between Zcd and Zci, 0.0005, 
however, is not distinguishable given the present accuracy 
of the energy scale. 

3.5. Other systematic errors 

The Point Spread Function (PSF) of the Suzaku X-ray 
telescopes has a width of ^ 2' (in half-power diameter), 
which is comparable to the size of the spectral regions 
examined here. This means that photons we detect in a 
certain spectral accumulation region are the sum of those 
coming from the corresponding sky region, Ci, and those 
from its surroundings, C2. We estimated this effect by 
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Fig. 4. Results of the redshift measurements using the He-hkc Fc-K lines in 52 cells. The left and right axes show the redshift and 
radial velocity of the ICM, respectively. Sec figure 2 for the definitions of the spectral regions. The error bars are la statistical. 
The optical redshift of z^i = 0.0104 and the range of the present calibrational error at the 90% confidence level, ±0.002 in redshift 
(or ±700 kms~^ in the line-of-sight velocity), are indicated with the horizontal dotted lines. 



the xissim ray-tracing simulator (Ishisaki et al. 2006), 
assuming a kT = 3 keV APEC model for the spectral dis- 
tribution and a double /3 model for the surface brightness 
distribution. According to Mohr et al. (1999), the two (3 
components arc assumed to have core radii of rd — 99 kpc 
and = 8.6 kpc, the same slope parameter (3 = 0.569, and 
a surface brightness ratio of /2//1 = 43 at the center. This 
simulation yielded Ci : C2 0.4 : 0.6 in cells lying at r < 4'. 5 
from the center, and Ci : C2 ^ 0.5 : 0.5 in cells surrounding 
them. Therefore, the results in figure 4 remains valid, as 
long as the spatial resolution of the analysis is regarded 
as 4' 50 kpc) instead of ~ 2'. 

As reported in Serlemitsos et al. (2006), the attitude of 
Suzaku sometimes drifts in both DETX/Y directions with 
a typical peak-to-peak amplitude about 0'.5 — 1', mainly 
synchronized with the spacecraft orbital period (96 min). 
This effect blurs the spatial photon distribution on the de- 
tector plane, because the size of the spectral accumulation 
regions is relatively small (~ 2') in our analysis. To exam- 
ine this issue, we performed Lorentzian fitting to central 
X-ray surface brightness profiles projected onto the DETX 
or DETY axis, and determined the X-ray centroid posi- 
tion of the cluster every 256 seconds, without relying on 
the attitude solution. As a result, the amplitudes of the 
drift were found to be small in this particular observation: 
ADETX - 37" and ADETY ~ 43". We produced a mock 
attitude file (that describes temporal change of the satel- 
lite's euler angles) assuming that the cluster centroid po- 
sition moves sinusoidally with a peak-to-peak amplitude 
of 40" in the DETX/Y directions and the 96-min period. 
By repeating the same ray-tracing simulation incorporat- 
ing this mock attitude file, the effect of the attitude drift 
on the photon counts (Ci or C2) was estimated to be 
< 10%, and typically Ci : C2 - 0.4 : 0.6 inside the XIS field 
of view. Therefore the effect of attitude drift is negligible 
compared to the original PSF effects. 

3.6. Analysis of larger-scale velocity structure 

In order to search the X-ray data for larger-scale veloc- 
ity gradients in the ICM, we compared the redshift in the 
following three regions; the central 18' x 18' square, and 
two 7' X 18' rectangular regions which are offset north- 
wards/southwards by about 12' (two dashed boxes in fig- 
ure 2). The spectral fits in the 5-10 keV band yield 



2 = 0.00993±0.00015(±0.0015),0.0119±0.0011 (±0.0015), 
and 0.0094± 0.0015 (±0.0015), for the central, northern 
offset, and southern offset regions, respectively. The 
68% statistical errors and (the 68% systematic errors) 
are quoted. The radial velocity is then Vr ~ — 144± 
44 (±470) kms-i, 441 ± 316 (±470) kms-\ and -306 ± 
450 (±470) kms~^ for the three regions, respectively. 
Thus, no significant velocity variation on 10-arcminutes 
scales is seen from the above analysis. Taking the sum of 
the statistical and systematic errors in quadrature, we es- 
timate the 90% upper limit on the velocity gradient among 
the above three regions to be |A?j| < 1400 kms"""^. 

4. Discussion 

From the careful analysis of the high-quality Fe-K line 
data obtained with the Suzaku XIS, we found no evidence 
of significant spatial gradients in the Fe-K line center ener- 
gies inside the XIS fields of view. The derived upper limits 
apparently exclude the large velocity difference previously 
claimed by Dupke & Bregman (2006). In §4.1, we further 
investigate this issue by analyzing the same regions as 
studied by Dupke & Bregman (2006). In §4.2, implica- 
tions of the present results are discussed, followed by §4.3 
where we constrain Doppler broadening of the Fe-K lines 
in search for turbulent motions of the ICM. 

4.1. Comparisons with previous results 

A large velocity difference by ~ 2500 kms~^ was re- 
ported in the core of the Centaurus cluster based on 
the previous ASCA and Chandra observations (Dupke 
& Bregman 2001; Dupke & Bregman 2005; Dupke & 
Bregman 2006). Specifically, by comparing the redshifts 
measured with the Chandra/ACIS-S3 spectra in various 
regions near the cluster center, Dupke & Bregman (2006) 
obtained (1) a maximum velocity gradient of Atj = 2900 ± 
700 kms~^ between a pair of rectangular regions, 2'. 3 x 
2'. 7 (PEAK-) and 1'.5 x 2'.5 (PEAK+), and (2) another 
large difference by Aw = 2400 ± 1000 kms~^ between two 
2'. 4 X 3' regions (R-2 and R+3). We show these four re- 
gions on our figure 2. 

Since the report by Dupke & Bregman (2006) appar- 
ently disagrees with our results, we examined the present 
Suzaku data for the above two cases, by accumulating 



No. ] 

spectra from the same sky regions on the three FI chips. 
(As seen in figure 2, the PEAK- and PEAK-I- regions are 
located near cells #35 and #f 8, and R-2 and R-l-3 are near 
cells #43 and #2f in om- analysis.) After co-adding the 
central and offset data, the total exposure time is 183 ks 
for PEAK- and R-2, and 188 ks for PEAK+ and R+3; 
these are about 5 times longer than the Chandra expo- 
sure. 

As to (1), the model fitting to the XIS spectra 
has given z_ = 0.0101 ± 0.0005 (±0.0015), and z+ = 
0.0079 ±0.0012 (±0.0015), resulting in a negative value 
for the redshift difference Az = z+ — z_ — —0.0022 ± 
0.0013 (±0.0022), or the velocity difference, Av = cAz = 
-660 ±390 (±660) kms"^ Because PEAK- and PEAK+ 
are close to each other, the intrachip gain difference be- 
tween them is considered to be smaller than our canonical 
value, ±470 kms~^. Nevertheless, to be conservative, we 
multiplied this by \/2, to obtain ±660 kms~^ as our sys- 
tematic error on the velocity difference. The above results 
then imply that the velocity difference between PEAK- 
and PEAK+ docs not significantly exceed the calibration 
uncertainty of ±660 kms~^. If we further sum the sys- 
tematic and statistical errors in quadrature, the 90% con- 
fidence range becomes —1890 kms~^ < Av < +570 kms~^, 
which is inconsistent with the Chandra measurement of 
2900 ±700 kms~^. The results arc summarized in table 2. 

As to (2), a similar analysis of the Suzaku XIS data 
has given the redshifts of the R-2 and R±3 regions as 
ZR_2 = 0.0107 ± 0.0006 (±0.0015) and zr+3 = 0.0089 ± 
0.0010 (±0.0015), respectively. Therefore, the redshift and 
velocity differences become Az = zr+s- zr_2 = — 0.0018 ± 
0.0012 (±0.0022) and Aw = -540 ± 360 (±660) kms"!, re- 
spectively (table 2). The allowed 90% range then becomes 
-1730 kms-i <Av< +650 kms-\ which again excludes 
the Chandra measurement of Ad between R-2 and R+3. 

Then, can we reconcile the apparently conflicting 
Chandra (plus ASCA) and Suzaku results? One possibil- 
ity is that the much wider PSF of the Suzaku XRT, com- 
pared to that of Chandra, has diluted the velocity gradi- 
ent. To examine the effect, we again executed the xissim 
simulation of the XIS spectra for PEAK- and PEAK+, 
simply assuming that the gas bulk motion exists only in 
those two regions and that they have the velocity differ- 
ence close to the 90% lower limit of the Chandra result, 
Av = +1700 kms~^; other regions in the XIS field of view 
are assumed to have no bulk motion, with z = 0.0104. 
The model fit to the two simulated spectra indicates that 
the Suzaku XIS should observe a velocity difference of 
Av = +690 kms~^, as a result of the telescope resolution 
which reduce the velocity difference by a factor of ~ 2.5. 
Therefore, the Chandra measurement may be brought to 
be marginally consistent with the present Suzaku upper 
limit. As to (2), the same simulation indicates that the 
90% lower-limit value of Av ~ +760 kms^^ between R-2 
and R+3 based on the Chandra report will be reduced to 
Av = +300 kms~^ in the XIS data after smeared by the 
PSF. This is consistent with the 90% upper limit set by 
Suzaku, +650 kms~^. 

In this way, the positive At; reported by Dupke & 
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Bregman (2006) for the two pairs of regions can be tenta- 
tively made consistent with the negative Aw we measured. 
However, this was possible by invoking rather extreme as- 
sumptions; namely, to employ the Chandra lower limit 
and the Suzaku upper limit on Av, and to assume that 
the ICM bulk motion is localized to the small regions of 
~ 25 kpc. 

A comparison of the present results with the re- 
ported ASCA measurement (Dupke & Bregman 2001) 
may give additional clues. The Suzaku XRT (Serlemitsos 
et al. 2006) is an improved version over its predecessor 
onboard ASCA. Similarly, the Suzaku XIS is an improved 
instrument based on the SIS experiment onboard ASCA. 
As compared with the XRT+SIS combination onboard 
ASCA, the XRT+XIS of Suzaku hence has a ~ 1.5 times 
better PSF, a ~ 4.5 times larger effective area (at 7 keV; 
two SIS sensors versus three XIS cameras). Furthermore, 
the energy scale of the Suzaku XIS has been better cal- 
ibrated (Koyama et al. 2006) than the ASCA SIS which 
was the first single-photon detection CCD used in cos- 
mic X-ray observations. Therefore, any Doppler shifts in 
the Fe-K line detectable by ASCA would be detected by 
Suzaku with a much higher significance (unless the expo- 
sure is too short). 

From these considerations, we conclude that it is rather 
difficult to reconcile the the present results with the 
Chandra (plus ASCA) measurements. Additional infor- 
mation will be obtained by analyzing a longer Chandra 
dataset on the Centaurus cluster. (Fabian et al. 2005). ^ 

4-2. Implications of the derived upper limits 

Utilizing the Fe-Ka line, we obtained a 90% upper limit 
on the velocity difference of |Af | < 1400 kms~^ between 
any pair of regions within the XIS field of view, sepa- 
rated typically by ~ 50 kpc or more. This gives a tighter 
constraint on the ICM bulk motion than was suggested 
by previous reports. Since the ion sound velocity in the 
ICM is s = (7P/p)i/2 = (7fcr/^mp)i/2 _ S80 kms^^ for 
a temperature of kT = 3 keV, our results imply that the 
line-of-sight velocity of the ICM bulk motion, if any, does 
not largely exceed the sound velocity. In other words, 
the ram pressure associated with the ICM bulk motion is 
at most comparable to the thermal ICM pressure, thus 
approximately validating the assumption of hydrostatic 
equilibrium in calculating the cluster mass. Of course, we 
can probe the gas motion only along the line of sight, and 
hence we need to assume 3-dimensional velocity structures 
in order to quantitatively utilize the present upper limits. 

As an elementary exercise, below we evaluate the up- 
per limit on rotational motion assuming for simplicity that 
the gas is rigidly rotating at a typical circular velocity of 
ar |Az;|/2, and derive its quantitative effect on clus- 
ter mass estimation. Assuming an approximate spherical 
symmetry, and including the centrifugal force, the bal- 
ance against the gravitational pull at a radius r on the 
rotational equatorial plane then becomes 

^ Some of the present authors actually analyzed this data set, but 
did not detect any significant velocity differences between the 
R-2 and R-l-3 regions identified by Dupke &; Bregman (2006). 
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Table 2. Comparisons of velocity differences measured with Suzaku and Chandra 



[Vol. , 



Region 


Suzaku/XIS 


Chandra/ACIS 




Av (kms^"'") 


Av (kms"^)f 


PEAK-, PEAK+ 


-660 ±390 (±660) 


2900 ±700 


R-2, R+3 


-540 ±360 (±660) 


2400 ±1000 



* The velocity difference derived from the XIS spectra. The 68% statistical errors and (the 68% systematic errors) are 
quoted. ^ The velocity difference and the la error derived from the Chandra ACIS spectra (Dupke & Bregman 2006). 



(1) 



j,2 Qj. J. 

where M(r) is the cluster mass interior to radius r, pgas 
is the gas density, Pgas = pgaskT{^mp)~^ is the ICM pres- 
sure, and / is the fraction of the ICM that is rotating (0 < 
/ < 1). The gas density profile of the Centaurus cluster 
roughly follows p. 



gas I 



within the central 10' region, 
assuming the double-/? model (Ikebe et al. 1999; Mohr et 
al. 1999). Together with the assumption Ur oc r, this yields 



(2) 



and hence the second term on the right-hand side of equa- 
tion 1 can be rewritten as 



fP. 



gasCT,, 



r dr 
Equation 1 then becomes 
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where 



''-^--(o^)( 



700 kms- 



kT 
3 kcV 



(5) 



is the ratio of the kinetic energy due to rotational motion 
to the thermal energy. Thus, the total mass should be 
higher than the hydrostatic mass by a factor of (1 + fPr), 
because of the presence of additional pressure support. 

At a typical radius of ~ 100 kpc, the rotational veloc- 
ity is ar < \Av\/2 = 700 kms-\ yielding {1 + fPr) < 2. 
Considering an inclination angle i of the rotation axis 
and the mean sini factor, 2/7r, the rotational velocity 
is ar < \Av\/2smi = 1100 kms"^ and (l + fPr) ^ 3.6. 
Therefore cluster mass estimation with the assumption 
of hydrostatic equilibrium is justified within a factor of 
~ 2 — 3.6. In order to put a stronger constraint on the 
cluster mass to compare with other techniques (partic- 
ularly gravitational lensing), a two- fold improvement on 
the energy-scale calibration is necessary 

4-3. Turbulent broadening of iron lines 

The velocity structure in the cluster can also be studied 
by analyzing the X-ray line widths for turbulent Doppler 
broadening of the X-ray emission lines. The iron ions 
provide the best tracers of this effect, because the broad- 
ening due to a turbulent velocity, cturb, on the order 
of ~ 300 kms~^, can exceed their thermal broadening 



(~ 2 eV) by several times (Inogamov & Sunyaev 2003). 
However, considering the energy resolution of the FI 
CCDs, 160 eV (FWHM) or 70 eV (Gaussian la), what 
can be done with the present data is at most to exam- 
ine the observed Fe-K line width for possible broadening 
beyond the instrumental energy resolution and its intrin- 
sic line width (ctq ~ 28 eV for kT = 3 keV). In order to 
carry out this search under the highest statistics, we again 
fitted the XIS0+XIS2+XIS3 spectrum from the central 
18' X 18' region (analyzed in §3.6) with the APEC plus 
three Gaussian model. Here we assumed that the gas has 
no bulk velocity, and fixed the centroid of the He-like Ka 
line to 6.677(1 + Zci)~^ keV. As a result, the Gaussian 
width has been obtained as a = 47 ± 2 eV. 

Since the degradation of the energy resolution with time 
is not currently implemented in the XIS response ma- 
trix (Koyama et al. 2006), the fitted Gaussian width is 
expressed by a^ = ^hist + '''o + ''"tmb ' where CTmst is the 
width of the Mn-Ka calibration line as mentioned be- 
fore, and the average of the three FI CCDs gives cTinst = 
40 ± 4 eV. Thus, the measured width is consistent with 
that expected from the instrumental and intrinsic widths, 
((Tf„,t + (7^)1/2 = 49 ± 2 eV, where the intr insic width ao 
includes also the thermal Doppler although it is negligi- 
bly small (~ 2 — 3 cV) in the energy range of interest. 
The turbulent broadening is then loosely constrained as 
CTturb < 32 eV - 1400 kms^^ (90% confidence). Although 
(To significantly increases below a temperature of 2 keV, 
it is within the range of 26 eV < ao < 30 eV as long as 
2 keV < kT < 4 kcV, and does not significantly affect the 
above discussion. The derived limit on crturb is consis- 
tent with a stronger constraint, < 800 kms~^, derived by 
an argument that strong turbulence would smear out, via 
diffusion, the observed abundance gradient in the central 
100 kpc of the cluster (Graham et al. 2006). 

Observations of resonance scattering of the X-ray emis- 
sion lines will provide another way to constrain the ICM 
turbulent motion. This is, however, beyond the scope of 
the present paper. 

5. Summary 

Based on Suzaku XIS observations of the Centaurus 
cluster, we investigated the bulk motion in the intraclus- 
ter medium to constrain its dynamical state. The high 
sensitivity and accurate calibration of the XIS sensors en- 
abled us to study positional dependence of the iron-K line 
energy on a 50 kpc scale, particularly over the central 
240 kpc region, as well as on a 140 kpc scale over a di- 
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ameter of 460 kpc. The resuhs show that there is no sig- 
nificant velocity gradient, within the present calibration 
uncertainty of ±660 kms~^, suggesting that the bulk ve- 
locity does not largely exceed the thermal velocity of the 
gas. The 90% upper limit on the velocity separation be- 
tween a pair of regions is \Av\ < 1400 kms~^, providing a 
tighter constraint than the previous observations by a fac- 
tor of ^ 2. The present result cannot be easily reconciled 
with the previous Chandra observations which reported 
the large velocity gradients in the two pairs of rectangu- 
lar regions near the cluster core. The mean rcdshift of the 
ICM is determined to be 0.0097, in agreement with the op- 
tical redshift of the cluster within the calibrational error. 
With the upper limit on the velocity difference, the pres- 
sure support due to hypothetical rotational motion may 
not alter the cluster mass estimate under the hydrostatic 
assumption by more than a factor of ^ 3. 

We are grateful to all the members of the Suzaku 
Science Working Group, especially T. Tsuru and H. 
Nakajima for comments on the XIS energy scale, Y. 
Ishisaki for useful suggestions on the ray-tracing simula- 
tions, and N. Y. Yamasaki for discussions. N. O. acknowl- 
edges support from the Special Postdoctoral Researchers 
Program of RIKEN. 
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